Temperature effects on the nuclear symmetry energy and symmetry free energy with 
an isospin and momentum dependent interaction 
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Within a self-consistent thermal model using an isospin and momentum dependent interaction 
(MDI) constrained by the isospin diffusion data in heavy-ion collisions, we investigate the tempera- 
ture dependence of the symmetry energy Esym{p,T) and symmetry free energy Fsym{p,T) for hot, 
isospin asymmetric nuclear matter. It is shown that the symmetry energy Esym{p,T) generally 
decreases with increasing temperature while the symmetry free energy Fsym{p,T) exhibits opposite 
temperature dependence. The decrement of the symmetry energy with temperature is essentially 
due to the decrement of the potential energy part of the symmetry energy with temperature. The 
difference between the symmetry energy and symmetry free energy is found to be quite small around 
the saturation density of nuclear matter. While at very low densities, they differ significantly from 
each other. In comparison with the experimental data of temperature dependent symmetry energy 
extracted from the isotopic scaling analysis of intermediate mass fragments (IMF's) in heavy-ion 
collisions, the resulting density and temperature dependent symmetry energy Esym{p,T) is then 
used to estimate the average freeze-out density of the IMF's. 

PACS numbers: 25.70.-z, 21.65.4-f, 21.30.Fc, 24.10.Pa 



I. INTRODUCTION 

The equation of state (EOS) of isospin asymmetric nu- 
clear matter, especially the nuclear symmetry energy, is 
essential in understanding not only many aspects of nu- 
clear physics, but also a number of important issues in as- 
trophysics [l|,[3,[l,[l,[l,i,0,B[l- Information about the 
symmetry energy at zero temperature is important for 
determining ground state properties of exotic nuclei and 
properties of cold neutron stars at /3-equilibrium, while 
the symmetry energy or symmetry free energy of hot 
neutron-rich matter is important for understanding the 
liquid-gas phase transition of asymmetric nuclear matter, 
the dynamical evolution of massive stars and the super- 
nova explosion mechanisms. Heavy-ion reactions induced 
by neutron-rich nuclei provide a unique means to investi- 
gate the symmetry energy [110, B]- In particular, recent 
analyses of the isospin diffusion data in heavy-ion reac- 
tions [13, [ll|, [ill have already put a stringent constraint 
on the symmetry energy of cold neutron-rich matter at 
sub-normal densities. On the other hand, the temper- 
ature dependence of the symmetry energy or symmetry 
free energy for hot neutron-rich matter has received so 
far little theoretical attention [l^, 0, [l^ . 

For finite nuclei at temperatures below about 3 MeV, 
the shell structure and pairing as well as vibrations of 
nuclear surfaces are important and the symmetry en- 
ergy was predicted to increase slightly with the increasing 
temperature [IB, [i3l ■ Interestingly, an increase by only 
about 8% in the symmetry energy in the range of T from 
to 1 MeV was found to affect appreciably the physics 
of stellar collapse, especially the neutralization processes 
[l6| . At higher temperatures, one expects the symme- 



try energy to decrease as the Pauli blocking becomes less 
important when the nucleon Fermi surfaces become more 
diffused at increasingly higher temperatures [l^, [13, [3l ■ 
Based on a simplified degenerate Fermi gas model at fi- 
nite temperatures, two of present authors [lH have re- 
cently studied the temperature dependence of the sym- 
metry energy and it was shown that the experimentally 
observed decrease of the nuclear symmetry energy with 
the increasing centrality or the excitation energy in iso- 
topic scaling analyses of heavy-ion reactions can be well 
understood analytically within the degenerate Fermi gas 
model. In particular, it was argued that the symmetry 
energy extracted from isotopic scaling analyses of heavy- 
ion reactions reflects the symmetry energy of hulk nu- 
clear matter for the emission source. Furthermore, it 
was found that the evolution of the symmetry energy ex- 
tracted from the isotopic scaling analysis is mainly due to 
the variation in the freeze-out density rather than tem- 
perature when the fragments are emitted in the reactions 
carried out imder different conditions. 

In the present work, within a self-consistent thermal 
model using an isospin and momentum dependent inter- 
action (MDI) constrained by the isospin diffusion data in 
heavy-ion collisions, we study systematically the temper- 
ature dependence of the nuclear matter symmetry en- 
ergy Esym{p,T) and symmetry free energy Fsym{p,T). 
It is shown that the nuclear matter symmetry energy 
Esym{p,T) generally decreases with increasing tempera- 
ture while the symmetry free energy Fgym {p, T) exhibits 
opposite temperature dependence. The decrement of the 
symmetry energy with temperature is essentially due to 
the decrement of the potential energy part of the symme- 
try energy with temperature. Furthermore, the difference 
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between the nuclear matter symmetry energy Egym {p, T) 
and symmetry free energy Fsy^ (p, T) is found to be quite 
small around nuclear saturation density. Using the re- 
sulting density and temperature dependent symmetry en- 
ergy Esym{p, T), we estimate the average freeze-out den- 
sity of the fragment emission source based on the mea- 
sured temperature dependent symmetry energy from the 
isotopic scaling analysis in heavy-ion collisions. 

The paper is organized as follows. In Section|IIl we in- 
troduce the isospin and momentum dependent MDI in- 
teraction and the detailed numerical method to obtain 
the EOS of the symmetric nuclear matter and pure neu- 
tron matter at finite temperatures. Results on the tem- 
perature dependence of the symmetry energy and sym- 
metry free energy are presented in Section IIIII In Sec- 
tion [iVl we discuss the experimental data of the isotopic 
scaling in heavy-ion collisions by means of the obtained 
density and temperature dependent symmetry energy. A 
summary is given in Section [V] 



II. HOT NUCLEAR MATTER EOS IN 
MOMENTUM DEPENDENT INTERACTION 

Our study is based on a self-consistent thermal model 
using a modified finite-range Gogny effective interaction, 
i.e., the isospin- and momentum-dependent MDI inter- 
action In the MDI interaction, the potential energy 
density V{p,T,5) of a thermal equilibrium asymmetric 
nuclear matter at total density p, temperature T and 
isospin asymmetry 5 is expressed as follows [ill; [3 1 



V{p,T,5) 
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In mean field approximation, Eq. ([1]) leads to the follow- 
ing single particle potential for a nucleon with momentum 
p and isospin r in the thermal equilibrium asymmetric 
nuclear matter, i.e., [HI. [Tsj 
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In the above t = 1/2 ( — 1/2) for neutrons (protons); 
a = 4/3; fri'^TP) is the phase space distribution func- 
tion at coordinate f and momentum p. The parameters 



Au{x),Ai{x),B,CT,TTCr-r and A have been assumed 
to be temperature independent and are obtained by fit- 
ting the momentum-dependence of U{p, T = 0, 5,p, r) to 
that predicted by the Gogny Hartree-Fock and/or the 
Brueckner-Hartree-Fock calculations, the zero tempera- 
ture saturation properties of symmetric nuclear matter 
and the symmetry energy of 31.6 MeV at normal nuclear 
matter density po = 0-16 fm""^ [l8| . The incompressibil- 
ity Kq of cold symmetric nuclear matter at saturation 
density po is set to be 211 MeV. The parameters Au{x) 
and Ai(x) depend on the x parameter according to 
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The different x values in the MDI interaction are in- 
troduced to vary the density dependence of the nuclear 
symmetry energy while keepin g ot her properties of the 
nuclear equation of state fixed |ll| and they can be ad- 
justed to mimic predictions on the density dependence of 
nuclear matter symmetry energy by microscopic and/or 
phenomenological many-body theories. The last two 
terms of Eq. ^ contain the momentum-dependence of 
the single-particle potential. The momentum dependence 
of the symmetry potential stems from the different in- 
teraction strength parameters Cr.-r and Cr,T for a nu- 
cleon of isospin r interacting, respectively, with unlike 
and like nucleons in the background fields. More specif- 
ically, we use Cr-T = -103.4 MeV and Cr,r = -11.7 
MeV. We note that the MDI interaction has been exten- 
sively used in the transport model for studying isospin ef- 
fects in intermediate e nergy h eav y -ion collisions induced 
by neutron-rich nuclei dl 111 [H, US iH, [H, [H, HI . In 
particular, the isospin diffusion data from NSCL/MSU 
have constrained the value of a; to be between and —1 
for nuclear matter densities less than about 1.2po [HI, [Hi, 
we will thus in the present work consider the two values 
of X = and x = —1. 

At zero temperature, /r(?^,p) = ■p-0(p/(T) — p) and 
the integral in Eqs. H]) and (O can be calculated analyt- 
ically [l8| . For an asymmetric nuclear matter at thermal 
equilibrium with a finite temperature T, the phase space 
distribution function becomes the Fermi distribution 
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where p.^ is the chemical potential and e{p,T,5,p,T) is 
the total single particle energy for a nucleon with isospin 
T and momentum p, which includes the kinetic energy 
and the single particle potential U{p, T, S,p^ r), i.e.. 
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with 
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where ri and T2 can be chosen as the same or different 
to mimic the last two terms of Eq. ^ . 

The chemical potential fir is therefore independent of 
the nucleon momentum p and can be determined from 



p^dp 



exp( 



(7) 



From Eq. (O, we can see that, at finite temperature, 
to obtain the chemical potential fXr requires knowing 
e(^p,T,S,p,T) (and thus Rr^^T^iP^P)) for all while 
from Eq. ^ knowing Rt-^^t2{p,p) needs further to know 
frif^TP) which again depends on the chemical potential 
Pr from Eq. (g]). Therefore. Eqs. (g]), (O, dH), and 
constitute closed sets of equations whose solution can be 
obtained by a self-consistency iteration, just as in the 
Hartree-Fock theory. 

Following the recipe used in Rcf. [2^, the self- 
consistency problem of Eqs. ([1]), (O, and ([7]) can 
be solved by the following iterative scheme. Firstly, we 
make an initial guess for RT^ r2{PiP) from the zero tem- 
perature condition, i.e.. 
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where p/(r) — h{3'K^ prY^"^ is the Fermi momentum. The 
right hand side of Eq. ([5]) is not related to the chemical 
potential p,r and thus the initial form of the single nu- 
cleon energy e^{p,T, S,p,t) can be obtained. Secondly, 
substitute e°(/9, T, 6,p, r) into Eq. ([7]) to obtain the initial 
chemical potential p^ for protons and neutrons. Then, 
use e'^ (p, T, S, p, r) and p^ to obtain new Rt^ (p, p) func- 
tion, namely, R\ .^^{p,p) from Eq. (|4]) and ([6]). This in 
turn gives the new single nucleon energy e^{p,T, S,p,t) 
from Eq. ([5]) and then the new chemical potential p]. can 
be obtained from Eq. ((T]). The cycle is repeated and a 
few iterations arc sufficient to achieve convergence for the 
chemical potential pr with enough accuracy. It should be 
mentioned that the neutron and proton chemical poten- 
tials are coupled with each other in asymmetric nuclear 
matter and thus the convergence condition must be sat- 
isfied simultaneously for neutrons and protons. 

With the self-consistency iteration, we can finally ob- 
tain the chemical potential pr and the single nucleon 
energy e{p,T,5,p,T) for an asymmetric nuclear matter 
at thermal equilibrium with a finite temperature T. The 
potential energy density V{p, T, 5) of the thermal equilib- 
rium asymmetric nuclear matter then can be calculated 
from Eq. ([T]) and the energy per nucleon E[p, T, 5) is then 
obtained as 



E{p,T,6) = - 
P 



/2 
d'p£^Jr{r,p} 



Furthermore, we can obtain the entropy per nucleon 
Sr (p, T, S) of the thermal equilibrium asymmetric nuclear 
matter as 

Sr(p,T,S) = / p^fn,- ln7iT- + (l-n^) ln(l-7iT-)](ip 

Ph'^ Jo 

(10) 

with the occupation probability 
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Finally, the free energy per nucleon F{p, T, 5) of the ther- 
mal equilibrium asymmetric nuclear matter can be ob- 
tained from the thermodynamic relation 



F(p, T, 5) = E{p, T, S) - Tj2Sr{p, T, S). 



(12) 
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FIG. 1: (Color online) Density dependence of the energy per 
nucleon E (left panel) and free energy per nucleon F (right 
panel) for symmetric nuclear matter at T = MeV, 5 MeV, 
10 MeV and 15 MeV with the MDI interaction. 

Using the MDI interaction, we can now calculate the 
energy per nucleon E(p, T, 5) and free energy per nucleon 
F(p, T, S) of nuclear matter at finite temperature from 
Eq. ^ and p^ . Shown in Fig. [l]is the density depen- 
dence of E(p,T,S) and F{p,T,S) for symmetric nuclear 
matter at T = MeV, 5 MeV, 10 MeV and 15 MeV using 
the MDI interaction with a; = and —1. For symmetric 
nuclear matter ((5 = 0), the parameter a; = would give 
the same results as the parameter x = — 1 as we have 
discussed above, and thus the curves shown in Fig. [T]are 
the same for x = and —1. From Fig. [1] one can see 
that the energy per nucleon E{p, T, S) increases with in- 
creasing temperature T while the free energy per nucleon 
F{p, T, 6) decreases with increment of T. The increment 
of the energy per nucleon E(p^ T, S) with the temperature 
reflects the thermal excitation of the nuclear matter due 
to the change of the phase-space distribution function 
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fr{r,p). With the increment of the temperature, more 
nucleons move to higher momentum states and thus lead 
to larger internal energy per nucleon. On the other hand, 
the decrement of the free energy per nucleon F{p,T,S) 
with T is mainly due to the increment of the entropy 
per nucleon with increasing temperature. This feature 
also implies that the increment of TS{p, T) with T is 
larger than the increment of E{p, T) with T. Further- 
more, the temperature effects are seen to be stronger at 
lower densities while they become much weaker at higher 
densities. At lower densities, the Fermi momentum ^^(t) 
is smaller and thus temperature effects on the energy per 
nucleon E{p, T, 6) are expected to be stronger. Mean- 
while, the entropy per nucleon becomes larger at lower 
densities where the particles become more free in phase 
space and thus leads to a smaller free energy per nucleon. 
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FIG. 2: (Color online) Same as Fig. [T]but for pure neutron 
matter using the MDI interaction with a; = (solid lines) and 
— 1 (dashed lines). 

Similarly, shown in Fig. [2] are the density dependence 
of the E{p, T, S) and F{p^ T, 6) for pure neutron matter 
at T = MeV, 5 McV, 10 McV and 15 MeV using the 
MDI interaction with a; = and —1. The temperature 
dependence of the E{p, T, S) and F{p, T, S) for pure neu- 
tron matter is seen to be similar to that of the symmetric 
nuclear matter as shown in Fig. [TJ However, the param- 
eters X = and —1 display different density dependence 
for the energy per nucleon E{p, T, d) and free energy per 
nucleon F{p,T,d), which just reflects that the parame- 
ters X = and —1 give different density dependence of 
the nuclear symmetry energy and symmetry free energy 
as will be discussed in the following. 



III. NUCLEAR SYMMETRY ENERGY AND 
SYMMETRY FREE ENERGY 



As in the case of zero temperature, phenomenologi 
cal and microscopic studies [l^ 13 
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FIG. 3: (Color online) E{p, T, 5) -E{p,T,5 = 0) as a function 
of 5"^ at temperature T = MeV (a), 5 MeV (b), 10 MeV (c) 
and 15 MeV (d) for three different baryon number densities 
p — 0.5p(), l.bpo and 2.5po using the MDI interaction with 
a; = 0. 



equation of state of hot neutron-rich matter at density 
p, temperature T, and an isospin asymmetry 6 can also 
be written as a parabolic function of S, i.e., 

E{p,T,d) ^ Eip,T,d = 0) + E,ymip,T)S^ + O{d^). (13) 

The temperature and density dependent symmetry en- 
ergy Esym{p, T) for hot neutron-rich matter can thus be 
extracted from E,y,r,{p, T) w E{p, T, (5 = 1) - E{p, T,S 
0). The symmetry energy Esymip,T) is the energy cost 
to convert all protons in symmetry matter to neutrons 
at the flxed temperature T and density p. In order to 
check the empirical parabolic law Eq. (fT^ for the MDI 
interaction, we show in Fig. ^E{p, T, S) - E{p, T,S ^0) 
as a function of 5^ at temperature T = MeV, 5 MeV, 
10 MeV and 15 MeV for three different baryon number 
densities p = 0.5po,l-5po and 2.5po using the MDI in- 
teraction with X = 0. The clear linear relation between 
E{p,T,S) - E{p,T,5 = 0) and 5^ shown in Fig. [3] indi- 
cates the validity of the empirical parabolic law Eq. 
for the hot neutron-rich matter. We note that the em- 
pirical parabolic law Eq. (fT3)) is also well satisfied for the 
parameter x = —\. 

Similarly, we can define the symmetry free energy 
Fsym{p,T) by the following parabolic approximation to 



the free energy per nucleon 
F(p,T,5)^F{p,T,8^Q)- 



-F,ym{p,T)S' + 0{S^). (14) 



The temperature and density dependent symmetry free 
energy Fsym [p, T) for hot neutron-rich matter can thus 
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FIG. 4: (Color online) Same as Fig. [3] but for the free energy 
per nucleon F{p,T, 5). 



F{p, T,S = 0) which is just the free energy cost to 
convert all protons in symmetry matter to neutrons at 
the fixed temperature T and density p. In order to 
check if the empirical parabolic law is also valid for the 
free energy per nucleon of hot neutron-rich matter, we 
show in Fig. 0] F{p,T,S) - F{p,T,S = 0) as a func- 
tion of (S^ at temperature T = MeV, 5 MeV, 10 MeV 
and 15 McV for three different baryon number densities 
p ~ 0.5pQ, 1.5po and 2.5po using the MDI interaction 
with X — 0. One can see from Fig. |4] that the parabolic 
law Eq. (fT4|) is also approximately satisfied though at 
low densities and high temperatures, the linear relation 
between F{p, T, S) - F{p, T,S = 0) and is violated 
slightly. For the parameter a; = — 1, we also obtained the 
similar conclusion. 

In Fig. [51 we show the density dependence of the sym- 
metry energy Esym{p,T) and the symmetry free energy 
Fsy,n{p,T) at T = MeV, 5 MeV, 10 MeV and 15 MeV 
using the MDI interaction with x = and —1. For dif- 
ferent choice of the parameter a: = and — 1 , Esym (p, T) 
and Fgym {p, T) display different density dependence with 
X = (—1) giving larger (smaller) values for the sym- 
metry energy and the symmetry free energy at lower 
densities while smaller (larger) ones at higher densities 
for a fixed temperature. Similar to the E(p, T, S) and 
F{p, T, S) as shown in Figs. [T] and [21 the temperature 
effects on the symmetry energy Esym{p, T) and the sym- 
metry free energy F^ym {p, T) are found to be stronger at 
lower densities while they become much weaker at higher 
densities. 

Interestingly, we can see from Fig. \5\ that the sym- 
metry energy Egy„i{p,T) and the symmetry free energy 
Fsym{PiT) exhibit opposite temperature dependence, 




FIG. 5: (Color online) Density dependence of the symmetry 
energy Esym{p,T) (left panel) and the symmetry free energy 
Fsym{p,T) (right panel) at T = MeV, 5 MeV, 10 MeV and 
15 MeV using the MDI interaction with a: = and —1. 



namely, with increasing temperature T, EsymipiT) de- 
creases while Fsymip, T) increases. This also means that 
Fsym{p,T) always has larger values than Esym{p,T) at 
fixed density and temperature since they are identical at 
zero temperature. At higher temperatures, one expects 
the symmetry energy Esym (p, T) to decrease as the Pauli 
blocking (a pure quantum effect) becomes less important 
when the nucleon Fermi surfaces become more diffused 
at increasingly higher temperatures [H, [13, [l3| . On the 
other hand, the symmetry free energy Fsym{p,T) is re- 
lated to the entropy per nucleon of the asymmetric nu- 
clear matter, which is not a pure quantum effect, and 
its increment with increasing temperature can be under- 
stood by the following expression 

Fsym (p, T) = Esym. {P, T) 

+T [Snip, T, 6 = 0) + Sp{p, T, d = 0)] 
-TS„{p,T,S = l). (15) 



The first term of the right hand side in Eq. p5)) is the 
symmetry energy Esym{p,T), which decreases with in- 
creasing temperature as discussed above. However, the 
total entropy per nucleon of the symmetric nuclear mat- 
ter is larger than that of the pure neutron matter and 
their difference becomes larger with increasing tempera- 
ture, which leads to a positive value for the difference 
between the last two terms of the right hand side in 
Eq. Therefore, Fsym{p,T) has larger values than 

Esy7n{p,T) at fixed density and temperature. Further- 
more, the increment of TS{p, T) with T is stronger than 
the increment of E{p, T) with T as mentioned above, and 
the combinational effects thus cause the symmetry free 
energy Fsym.{p,T) increase with increasing temperature. 

Within the present self-consistent thermal model, be- 
cause the single particle potential is momentum depen- 
dent with the MDI interaction, the potential part of the 
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symmetry energy is expected to be temperature depen- 
dent. It is thus interesting to study how the potential and 
kinetic parts of the symmetry energy Esym (p, T) may 
vary respectively with temperature. However, for the 
symmetry free energy Fgym {p, T) , one cannot separate 
its potential and kinetic parts since Fgym {p, T) depends 
on the entropy that is determined by the phase space 
distribution function. Fig. [S] displays the temperature 
dependence of the symmetry energy Esym (p, T) as well 
as its potential and kinetic energy parts using the MDI 
interaction with x = at p = poi O.Spo: E^d O.lpo- With 
the parameter a; = — 1, the same conclusion is obtained. 
It is seen that both the symmetry energy Egym {p, T) and 
its potential energy part decrease with increasing temper- 
ature at all three densities considered. While the kinetic 
energy part of the Esym {p, T) increases slightly with in- 
creasing temperature for p ^ po and 0.5po a-nd decreases 
for p = O.lpo- These features are uniquely determined 
by the momentum dependence in the MDI interaction 
within the present self-consistent thermal model. The 
decrement of the kinetic energy part of the symmetry en- 
ergy with temperature at very low densities is consistent 
with predictions of the Fermi gas model at high temper- 
atures and/or very low densities [l^, H^. In the study of 
Ref. [iBl using the simplified degenerate Fermi gas model 
the potential part of the symmetry energy was assumed 
to be temperature independent for simplicity. The de- 
crease of the symmetry energy observed there is thus 
completely due to the decrease in the kinetic contribu- 
tion. However, we note that the temperature dependence 
of the total symmetry energy Esym{p,T) is consistent 
with each other for the two models. This is due to the 
fact that the phase space distribution function will vary 
self-consistently according to if the single particle poten- 
tial is or not momentum dependent. From the present 
self-consistent thermal model with momentum dependent 
MDI interaction, our results indicate that the decreasing 
symmetry energy with increasing temperature is essen- 
tially due to the decrement of its potential contribution. 



IV. ISOTOPIC SCALING IN HEAVY-ION 
COLLISIONS 

It has been observed experimentally and theoretically 
in many types of reactions that the ratio R2i{N,Z) of 
yields of a fragment with proton number Z and neu- 
tron number N from two reactions reaching about the 
same temperature T satisfies an exponential relationship 

EoifiV, zioc exp(aAf) [13, M, M, M, M, M, HI, H, M, 

l36l . 1371 . |38|. Particularly, in several statistical and dy- 
namical models under some assumptions [s^ . Issl . [s^ , it 

has been shown that the scaling coefficient a is related 
to the symmetry energy Cgym {p, T) via 
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FIG. 6: (Color online) Temperature dependence of the sym- 
metry energy Esym.{p,T) as well as its potential energy part 
and kinetic energy part using MDI interaction with s = at 
p = po (a), 0.5po (b), and O.lpo (c). 



where A[{Z/A)^] = [Zxj A^f -[Z2I is the difference 
between the {Zj A)^ values of the two fragmenting sources 
created in the two reactions. 

As mentioned in Ref. [l^, however, because of the dif- 
ferent assumptions used in the various derivations, the 
validity of Eq. (fTB|) is still disputable as to whether 



^Csyra{p. T) 



A [{Z/Af 



(16) 



and when the Csym, is actually the symmetry energy or 
the symmetry free energy. Moreover, the physical inter- 
pretation of the CsyrniPiT) is also not clear, sometimes 
even contradictory, in the literature. The main issue is 
whether the Csym measures the symmetry energy of the 
fragmenting source or that of the fragments formed at 
freeze-out. This ambiguity is also due to the fact that 
the derivation of Eq. pB]) is not unique. In particular, 
within the grand canonical statistical model for multi- 
fragmentation [3^ . [3^ the Csym refers to the symmetry 
energy of primary fragments. While within the sequen- 
tial Weisskopf model in the grand canonical limit [sJl 
it refers to the symmetry energy of the emission source. 
Following the arguments in Ref. (isj . we assume in the 
present work that the Csym reflects the symmetry energy 
of bulk nuclear matter for the emission source. 

In Fig. [3 we show the symmetry energy Esym{p,T) 
and symmetry free energy Esym {p, T) function of 

temperature using MDI interaction with x = and — 1 
at different densities from O.lpo to po- The temperature 
dependence of the symmetry e nerg y Egymip, T) is seen to 
be very similar to that in Ref. ^M, where a simplified de- 
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FIG. 7: (Color online) Temperature dependence of the sym- 
metry energy (solid lines) and symmetry free energy (dashed 
lines) using MDI interaction with 2; = (left panel) and 
— 1 (right panel) at different densities from O.lpo to po. 
The experimental data from Texas A&M University (solid 
squares) and the INDRA-ALADIN collaboration at GSI 
(open squares) are included for comparison. 



generate Fermi gas model at finite temperatures has been 
used. The symmetry energy does not change much with 
the temperature at a given density, especially around the 
saturation density po- Furthermore, it is seen that while 
the symmetry energy Esym{p,T) deceases slightly with 
the increasing temperature at a given density, the sym- 
metry free energy Fsym{p,T) increases instead. Around 
the saturation density po, it is found that the difference 
between the symmetry energy Egym (p, T) and the sym- 
metry free energy Fsymip, T) is quite small, i.e., only sev- 
eral percents, though at higher temperature, compared 
with their values at T = MeV. This feature confirms the 
assumption on identifying Cgym (p, T) to E,^yrn.(p, T) at 
lower temperatures and not so low densities [23, [30] ■ 
At low densities, on the other hand, the symmetry free 
energy Fsymip,T) exhibits a stronger temperature de- 
pendence and it is significantly larger than the symmetry 
energy Esym{p,T) at moderate and high temperatures. 
This is due to the fact that the entropy contribution to 
the symmetry free energy Fsymip^T) becomes stronger 
at low densities as mentioned in Sec. HI] and Sec. IIIII It 
should be noted that, at low densities the entropy may be 
affected strongly by the clustering effects [H, Ull which 
are not included in the present work. 

Experimentally, the temperature T and scaling coeffi- 
cient a (thus the Csym) of the fragment emission source 
can be directly measured while the determination of the 
density of emission source usually depends on the model 
used. Also included in Fig. [7] are the experimental data of 
the measured temperature dependent symmetry energy 
from Texas A&M University (TAMU) (solid squares) [30| 
and the INDRA-ALADIN collaboration at GSI (open 



squares) [3l|, d^]. From Fig. [71 it is seen clearly that 
the experimentally observed evolution of the symmetry 
energy is mainly due to the change in density rather than 
temperature, as shown in Ref. [l5|. Meanwhile, we can 
estimate from Fig. [7]the average freeze-out density of the 
fragment emission source from the measured temperature 
dependent symmetry energy based on the isotopic scaling 
analysis in heavy-ion collisions. In particular, using the 
symmetry energy E^ym (p, T) from the MDI interaction 
with a; = 0, we find the average freeze-out density of the 
fragment emission source p/ is between about 0.41po and 
0.52po for TAMU data while about 0.42po and 0.75po 
for INDRA-ALADIN collaboration data. On the other 
hand, using the symmetry energy Esymip,T) from the 
MDI interaction with a: = —1, the p/ is found to be 
between about 0.57po and 0.68po for TAMU data while 
about 0.58po and 0.84po for INDRA-ALADIN collabora- 
tion data. It is interesting to see that the extracted values 
of pf from the MDI interaction with a; = is very similar 
to those extracted in Ref. [13] using different models. 

Furthermore, if the symmetry free energy Fsym{p,T) 
from the MDI interaction with x = is used to estimate 
the pf, we find the p/ is between about 0.36po and 0.49po 
for TAMU data and about 0.33po and 0.72po for INDRA- 
ALADIN collaboration data. While if the symmetry free 
energy Fsym{p, T) from the MDI interaction with x = —1 
is used, the p/ is between about 0.52po and 0.66po for 
TAMU data and about O.Slpo and 0.83po for INDRA- 
ALADIN collaboration data. Therefore, the extracted pj 
values are not sensitive to if the measured Csym{p,T) is 
the symmetry energy or the symmetry free energy. How- 
ever, the extracted p/ values are indeed sensitive to the 
X parameter used in the MDI interaction, namely, the 
density dependence of the symmetry energy. We note 
that the zero-temperature symmetry energy for the MDI 
interaction with a; = and —1 can be parameterized, 
respectively, as 31.6(p/po)°-^^ MeV and 31.6(p/po)^ °^ 
MeV Therefore, the isotopic scaling in heavy- ion 

collisions provides a potential good probe for the density 
dependence of the nuclear matter symmetry energy once 
the average density of the emission source has been de- 
termined in the isotopic scaling measurement, as pointed 
out in Ref. [H. 



V. SUMMARY 

Within a self-consistent thermal model using the 
isospin and momentum dependent MDI interaction with 
a; = and —1 constrained by the isospin diffusion data 
in heavy- ion collisions, we have investigated the temper- 
ature dependence of the nuclear matter symmetry en- 
ergy Esym{p,T) and symmetry free energy Fsym{p,T). 
It is shown that the nuclear matter symmetry energy 
Esym{p,T) generally decreases with increasing tempera- 
ture while the symmetry free energy FsymiPi T) exhibits 
opposite temperature dependence. The decrement of the 
symmetry energy with temperature is essentially due to 
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the decrement of the potential energy part of the symme- 
try energy with temperature. The temperature effects on 
the nuclear matter symmetry energy and symmetry free 
energy are found to be stronger at lower densities while 
become much weaker at higher densities. Furthermore, 
the difference between the nuclear matter symmetry en- 
ergy Esymip, T) and symmetry free energy Fsymip, T) is 
found to be quite small around nuclear saturation den- 
sity, although significantly large at very low densities. 

Comparing the theoretical density and tempera- 
ture dependent symmetry energy Egy^ [p, T) with the 
Csym{p,T) parameter extracted from the isotopic scal- 
ing data from TAMU and the INDRA-ALADIN collab- 
oration at GSI, we found that the experimentally ob- 
served evolution of the symmetry energy is mainly due 
to the change in density rather than temperature, as 
shown in the previous work p^ . Meanwhile, we have 
estimated the average freeze-out density of the fragment 
emission source formed in these reactions by comparing 
the calculated Esym{p,T) or Fsym{p,T) with the mea- 
sured Csym (p, T) . Our results indicate that the extracted 
average freeze-out densities are not sensitive to whether 
the experimentally measured Cgym (p, T) parameter is the 
symmetry energy Esym {p, T) or the symmetry free en- 



ergy Fsym{p,T) in the temperature and density ranges 
reached in the TAMU and INDRA/GSI experiments. 
They are, however, sensitive to the x parameter used 
in the MDI interaction, namely, the density dependence 
of the symmetry energy. Therefore the isotopic scaling 
in heavy-ion collisions provides a potentially good probe 
for the density dependence of the nuclear matter symme- 
try energy provided the average density of the emission 
source can be determined simultaneously in the isotopic 
scaling measurements. 
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